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^ ' Abstract 

in ■ 

■rj" ■ We describe a quantum algorithm to prepare an arbitrary pure state of a register 

, of a quantum computer with fidelity arbitrarily close to 1. Our algorithm is based on 

Q I Grover's quantum search algorithm. For sequences of states with suitably bounded 

' amplitudes, the algorithm requires resources that are polynomial in the number 

of qubits. Such sequences of states occur naturally in the problem of encoding a 
classical probability distribution in a quantum register. 



g '■ 1 Introduction 



In many applications of quantum computers, a quantum register, composed of a fixed 
^ , number of qubits, is initially prepared in some simple standard state. This initial prepa- 
^ I ration step is followed by a sequence of quantum gate operations and measurements. 

There are applications of quantum computers, however, notably the task of simulating 
the dynamics of a physical system 13 El > that may require the initialization of a quan- 
tum register in a more general state, corresponding to the initial physical state of the 
simulated system. This leads naturally to the question of what quantum states can be 
efficiently prepared on a quantum register. 

The memory of a classical computer can be easily put into any state by writing an arbitrary 
bit string into it. The situation for quantum computers is very different. The Hilbert 
space associated with a quantum register composed of as few as 100 quantum bits (qubits) 
is so large that it is impossible to give a classical description of a generic state vector, i.e., 
to list the 2^°° complex coefficients defining it. In this sense it can be said that arbitrary 
pure states cannot be prepared jlj. 

It is nevertheless possible to formulate the problem of arbitrary state preparation for a 
register of qubits in a meaningful way. This is achieved by starting from the assumption 
that the state is initially defined by a set of quantum oracles. By assuming that the state 
is given in this form, we shift the focus from the problem of describing the state to the 
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problem of the computational resources needed to actually prepare it. In other words, 
we address the computational complexity rather than the algorithmic complexity of the 
state. For the purpose of quantifying these computational resources, we simply count the 
number of oracle calls. We are thereby deliberately ignoring the internal structure of the 
oracles and the computational resources needed in each oracle call. The algorithm we 
describe here is applicable to any set of oracles, i.e., to any state. We will show that it 
is efficient for a large class of states of particular interest for the simulation of physical 
systems. 

Let iV be a positive integer. We will describe a quantum algorithm for preparing a 
[log2 iVl-qubit quantum register in an approximation to the state 

N-l 

\^) = J2^/^e^^'^(^^\x) (1) 

x=0 

for arbitrary probabilities p{x) and arbitrary phases Here and throughout the paper, 
|0), . . . denote computational basis states. More precisely, given any small positive 
numbers A and u, our algorithm prepares the quantum register in a state |^) such that, 
with probability greater than 1 — z/, the fidelity obeys the bound 

1(^1^)1 >1- A. (2) 

To define the algorithm and to assess its efficiency for large N, we need to specify in 
which form the coefficients p{x) and are given. We assume that we are given classical 
algorithms to compute the functions p{x) and for any x. These classical algorithms 
are used to construct a set of quantum oracles. We will quantify the resources needed by 
our state preparation algorithm in terms of (i) the number of oracle calls, (ii) the number 
of additional gate operations, and (iii) the number of auxiliary qubits needed in addition 
to the [log2 A^l register qubits. 

To analyze the asymptotic, large N, behavior of our algorithm, we consider a sequence of 
probability functions pat : {0, . . . , A^ — 1} — > [0, 1], J^xP'^i-'^) ~ ^ sequence of phase 

functions (pN '■ {0, . . . ,N — 1} — > [0, 1], where N = 1,2,.... For any A^, the algorithm 
prepares the quantum register in a state such that, with probability greater than 
1 — z/, the fidelity obeys the bound 0, where in the definition (H)) of the functions p 
and are replaced hj pn and (pN, respectively. Under the assumption that there exists a 
real number rj, < r] < 1, such that 

Pn{x) < — ^ for all and x , (3) 

we show that the resources needed by our state preparation algorithm are polynomial in 
the number of qubits, log2 A^, and the inverse parameters ri~^, and u^^. 

An obvious example of a sequence of functions that do not satisfy the bound ^ and for 
which the resources required for state preparation scale exponentially with the number 
of qubits is given by Pn{x) = 6xy for some integer y = y{N). In this case, it follows 
from the optimality of Grover's algorithm 15] that the number of oracle calls needed is 
proportional to \/N. 
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Sequences that do satisfy the bound Q arise naturally in the problem of encoding a 
bounded probability density function / : [0, 1] [0, /max] in a state of the form 

N-l 

\^f)=Ar-'J2,/fWN)\x), (4) 

x=0 

where A/" is a normalization factor. Grover and Rudolph have given an efficient algorithm 
for this problem if the function / is efficiently integrable [7j. Essentially the same al- 
gorithm was found independently by Kaye and Mosca jS], who also mention that phase 
factors can be introduced using the methods discussed in Ref . [9] . Recently, Rudolph (TUI 
has found a simple nondeterministic state-preparation algorithm that is efficient for all 
sequences satisfying the bound Q- 

For general sequences of states satisfying the bound Q, a given value of the fidelity bound 
A, and assuming polynomial resources for the oracles, our algorithm is exponentially more 
efficient than the algorithm proposed by Ventura and Martinez ^Tj and later related 
proposals f2lEl; for which the resources needed grow like iVlog2 A^. The use of Grover's 
algorithm for state preparation has been suggested by Zeng and Kuang for a special class 
of coherent states in the context of ion-trap quantum computers ^ general analysis 
of the state preparation problem in the context of adiabatic quantum computation was 
given by Aharonov and Ta-Shma [T^ . 

This paper is organized as follows. In Sec. Elwe give a full description of our algorithm. A 
detailed derivation is deferred to Sec. 121 The algorithm depends on a number of parameters 
that can be freely chosen. In Sec. 0] we consider a particular choice for these parameters 
and show that it guarantees the fidelity bound (j21). We use the same choice of parameters 
in Sec. El to derive worst case bounds on the time and the memory resources required by 
the algorithm. In Sec. IHlwe conclude with a brief summary. 



2 Algorithm 

Our algorithm consists of two main stages. In the first stage, the algorithm prepares the 
register in an approximation to the state 

N-l 

l^p) = 5^ VpMk) , (5) 

x=Q 

which differs from only in the phases More precisely, let e be the largest small 

parameter such that 

e < Ar//3 (6) 
and 1/e is an integer. The first stage of the algorithm prepares the register in a state 

N-l 

l^p) = , (7) 

x=0 

such that, with probability greater than 1 — u, we have that 

\{^,\^,)\>1-X. (8) 
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We will describe the details of the first stage below. 

The second stage of the algorithm adds the phases to the state |\l/p) resulting from 
the first stage. This can be done in a straightforward way as follows. We start by choosing 
a small parameter e' such that 1/e' is a positive integer. We then define a list of unitary 
operations, Ui, . . . , Ui/^', on our quantum register by 



[_ \x} otherwise. 

The operators Uk are conditional phase shifts that can be realized as quantum gate se- 
quences using the classical algorithm for computing the function (j){x) [9j. If we apply the 
operators Uk sequentially to the result of the first stage, we obtain 
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'\x) if > {k- l)e' , 



^) = U^U, ■ ■ ■ f/i/e'l^p) = V e2-^^(^')|x) , (10) 




where the function 0(x) satisfies the inequality 

|0(x)-0(x)| <e72 (11) 

for all X. It can be shown (see section 14. 4p that together with Eq. (jSI) this implies the 
bound 

1(^1^)1 > 1 -A- A' , (12) 

where A' = e'^/8. Notice the slight abuse of notation identifying the parameter A in the 
inequality with the sum A + A' in the inequality (fT^ . 

We now proceed to a more detailed description of the first stage of the algorithm. From 
now on we assume that N is an integer power of 2. This can always be achieved by 
padding the function p{x) with zeros. Given our choice of the parameter e, Eq. (jH)), we 
define a list of oracles, Oi, . . . , Oi/^, by 



oJx) = i^ if VpM>(l-eA;)/v^, .^3. 
\ otherwise. 

We extend this definition beyond the domain of the function p by setting Ok{x) = for 
X > N. Using the classical algorithm to compute p{x), one can construct quantum circuits 
implementing the unitary oracles 

6fc|a:) = (-l)°^-(^)|a;). (14) 

These circuits are efficient if the classical algorithm is efficient. The list of oracles 
defines a new function p'{x) via 

^p'{x) = (1 - ek)/^^ if Ok-i{x) = and Ok{x) = 1 , (15) 

where Oo{x) = by convention. The situation is illustrated in Fig. d where the x values 
have been permuted for clarity. Knowledge of this permutation is not required for our 
algorithm. 



4 




Figure 1: (Color online) The solid line shows an example for a function ^/p{x), sketched 
versus a permutation, n(x), chosen so that p(Il~^{x)) < p(Il~^{y)) if x > y. The dotted 
and dashed lines show the corresponding functions \/p'{x) and \/p"{x) [see Eqs. f|TH|) 
and (f^H|) ]. which in this representation look like decreasing step functions. In this example, 
1/e = 7, T = 4, and it is assumed that < due to counting errors. Notice that the 
function p{x) is not required to be decreasing. 



The essence of the first stage of the algorithm consists in using a number of Grover 
iterations based on the oracles dk to prepare the register in an approximation to the state 

N-l 

\m^,)=M'-'Y.^^)\x) , (16) 

where the normalization factor M' reflects the fact that p'{x) may not be normalized. To 
find the number of required Grover iterations for each oracle 6^, we need an estimate of 
the number of solutions, ra^, for each oracle, defined by 

N-l 

rik = ^Ok{x) . (17) 

x=0 

This estimate can be obtained from running the quantum counting algorithm [H] for each 
oracle dk- We denote the estimates obtained in this way by hk- The accuracy of the 
estimate relative to Uk can be characterized by two real parameters, < u < 1 and 
< rjc < 1/2, in such a way that, as a result of quantum counting, with probability 
greater than 1 — z/ we have 

\nk - nk\ < VcN for k = 1, . . . ,1/e . (18) 

For each oracle, 6k, the resources needed to achieve the counting accuracy specified by 
rjc and u depend on the actual number of solutions Uk- This dependence is important 
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for optimizing the performance of our algorithm. In this paper, however, we present a 
simpler analysis assuming worst case conditions for each oracle Ofc. For this analysis we 
use a specific choice of ?7c, which is given by Eq. 

The analysis of the algorithm is simplified if we concentrate on a subset of oracles, 

Ok = Of^ for k = l,...,T, (19) 

where T and the indices f\, . . . , fx are determined by the construction below. We intro- 
duce a new parameter rjg (see Eq. (j^^ below) such that rjc < rjg < 1/2. The index /i is 
defined to be the smallest integer such that 

^/i > Vg^ fij < TjgN for j < fi , (20) 

and, for k >2, the index fk is the smallest integer such that 

f^fk-i < and fk-i < fk< 1/e • (21) 

The number T is the largest value of k for which these inequalities can be satisfied. The 
effect of Eq. fl20|) is to neglect narrow peaks (corresponding to small values of 11 (x) in 
Fig.P). Equation ()21|1 makes sure that the numbers n/^. form an increasing sequence even 
if, due to counting errors, the estimates hk do not (see Fig. [T)). 

For k = 1, . . . ,T, define Nk = 

Af-l 

Nk = J20k{x), (22) 

and 

h = e{h+i-fk)/VvN, (23) 

where we define /t+i = 1/^- For every oracle Ok, the value of Nk is an approximation to 
the number of solutions, Nk, satisfying the bound 

\Nk-Nk\ = \nf,-nf^\<r],N. (24) 
In what follows it will be convenient to introduce the notation 

k 

Bs,k = 5^ 5, . (25) 

j=s 

The oracles Ok define a new function p"{x) via 

= Bk,T if Ok-i{x) = and Ok{x) = 1 , (26) 

where Oo{x) = by convention. The function ^yp"(Il~^(x)) is a decreasing step function, 
with step sizes Si, . . . ,St which are multiples of e/ y/rjN. The widths of the steps are given 
by the numbers Nk which are determined by the oracles (see Fig. [T)) . 
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The algorithm can now be completely described as follows. Choose a suitable (small) 
number, a, of auxiliary qubits (see Eq. (jU^ below), and define L = log2 + a. For 
k = 1, . . . ,T, find the quantities 



k-l „ fe-l 



(^iVA) +iV,(l-5^(iV,-iV,_i)i?,Vi 



and 



al = ^ . 2=1^ , (27) 

Nk{2-N-Nk) 

~fin^ and 7r = , (28) 

f 2Nk\ 

Cjh = arccos 1 , (29) 

(30) 
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For /c = 1, . . . , T, define the Grover operator 

GiPuM) = ((2|*°)(*°| , (31) 

where 

1^0) = (2'^Ar)-i/2 J2 \x) , (32) 

/ is the L-qubit identity operator, Ok\x) = (— 1)*^'='-^'' and where the domain of the 
oracles Ok is extended to the range < x < 2°-N by setting Ok{x) = ii x > N. 

Prepare a register of L qubits in the state I^P^), then apply the Grover operators succes- 
sively to create the state 

|M/^) = G(OT,tT)---G(Oi,ti)|vl/0). (33) 

Now measure the a auxiliary qubits in the computational basis. If all a outcomes are 
this stage of the algorithm is successfully prepares the desired state Eq. ((7j). 

If one of the measurements of the auxiliary qubits returns 1, this stage of the algorithm 
has failed, and one has to start again by preparing the register in the state 1^^°) as in 
Eq. (1221) • Assuming the choice of parameters in Eq. the probability, p^i, that the 
algorithm fails in this way satisfies the bound pfaii < lOA (see Eq. p32|l ). 

Before we provide detailed proofs of the above claims it is helpful to give a hint of how 
this stage of the algorithm achieves its goal. The algorithm aims at constructing the 
function \/p{x) which is close to the function \/p"{x) defined in Eq. (jSH))- The sequence 
of Grover operators in Eq. creates a step function that is close to ^y p"{Il~^{x)) (see 
Fig. H}. In particular each operator G{Ok,tk) in Eq. (jHH)) creates a step with the correct 
width Nk and a height hk which is close to the target height 6k. Due to a remarkable 
property of Grover 's algorithm (THj, once the features hi, ... , hk-i have been developed 
they are not distorted by G{Ok, tk) which develops hk. In this way the algorithm proceeds 
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building feature after feature until it constructs all T of them. At the end, because of the 
inherent errors, the auxiliary qubits end up having small amplitudes for nonzero values. 
Measuring them projects the auxiliary qubits onto the zero values with a probability that 
can be made arbitrarily close to 1. This also slightly changes the features hk due to 
renormalization of the state after the measurement, which we take into account when we 
estimate the overall loss of fidelity. 



3 Derivation of the algorithm 

In section El we have already explained the second stage of our algorithm. Here we present 
a detailed explanation of the first stage. This section is organized as follows. In subsec- 
tion we review some properties of the Grover operator introducing our notation as we 
go along. In subsection 13.21 we introduce a convenient mathematical form for analyzing 
intermediate quantum states visited by the algorithm. And finally, in subsection 13.31 we 
derive the values PU)) of the times tk used in our algorithm. 



3.1 Preservation of features by the Grover operator 

We will be using the following result [15]. Consider an oracle O, which accepts r values 
of X (out of the total of 2°'N, i.e., 0{x) = r). We shall call such values of x good, 

as opposed to bad values of x that are rejected by the oracle. Using different notation for 
the coefficients of good and bad states, we have that after t Grover iterations an arbitrary 
quantum state 

1^'") = E + E (34) 

good X bad x 

is transformed into 

|vl/fi-) = G{0,t)\¥^) = E 9T\^) + E ^^1^)- (35) 

good X bad x 

Let g™^ and 6'" be the averages of the initial amplitudes of the good and the bad states 
respectively: 



-im 



r ^ ' 2^N-r 

good X bad x 



and similarly for the final amplitudes 



r- = ^E^^' ^''^ = ^^E^^' (37) 

good X bad x 



Let us also define 

Agf = gf - , = U^' - U^' . (38) 
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In other words, A^f™ and A6™ define the features of the initial amplitude functions ^f™' 
and 6™^ relative to their averages and 6™'. Biham et. al. have shown that the change 
of the amplitudes is essentially determined by the change of the averages: 

= fefi-^ + (-l)*A6^" , (39) 
where the averages g^^ and are given as follows. Define 

2r 

uj = arccos 1 — 



a = W Ift'-^'P + Ifli^P — 



^ini y - r 



arctan(£-J— — - ) . (40) 



The averages are given by 



— r 

g^"" = \ a sm{ujt + cp) , 

fefi" = acos(cjt + 0). (41) 
We shall also use the separations, A"^' and A*^°, of the averages 

^ini _ -ini _ ^ini ^fin _ -fin _ ^fin _ ^42) 

Directly from the definition we obtain 

A^'^ = ^^^sin(cjt-0, (43) 

where the phase can be found from A™', 



e = -arcsin —./-—). (44) 



a V 2"A^ 

In our applications we will only need the case when the initial amplitude of the bad states 
is fiat, i.e. A6^°' = 0. In this case, the amplitude of the bad states always remains fiat 

6f = 6^" , (45) 

This fact and the fact that the initial features, A^f™^ of the amplitude of the good states 
are preserved, 

^fin _ -fin ^ ^^ini ^ (4q) 

see fj39|l . is crucial for understanding the rest of this paper. 
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x„ X x„ X 



Figure 2: (Color online) Illustration of the action of the Grover operator G{0,t) on an 
arbitrary amplitude function ^/p{x)■ The light area indicates values of x that are accepted 
by O ("good values"), and the dark area indicates values of a; that are rejected by O ("bad 
values" ) . 

3.2 Intermediate states 

The preparation stage of our algorithm summarized in Eq. ()33|) gives rise to a sequence 
of states defined by 

\^'') = G{Ok,tk)\^''') . (47) 
Let Ofe be the set of solutions to the oracle Ok- 

Ok = {x:Okix) = l}. (48) 

For k = 2, . . . ,T, these states can be written in the form 

|vI/^) = At ^ + + + (49) 

where 

k 

4 = B'' + J2hs- (50) 

s=j 

Since is real and positive, the value of B'^ can be determined from the normalization 
condition = 1. 

The action of the algorithm can be visualized as shown in Fig. El which shows the result 
of the first three iterations. The integers [defined in Eq. ((221)] are the number of good 
values of x according to oracle Ok- We see that each operation G{Ok, tk) prepares a feature 
of height hk = hk{tk, Ni, . . . , Nk)- It follows from the conclusions of section IXD that once 
such a feature is developed, its height remains constant throughout the computation. 

3.3 Values of tk 

In this subsection we show how the times tk are related to the corresponding features hk- 
The normalization condition (\l>'^|vl>'^) = 1 reads 

NMiY + - N,)iA',y + . . . + (AT, - Nk.i){AlY + (2'^iV - Nk)iB>'Y = 1. (51) 
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q(x) 



^3 + 

53 4 



hsits,N,,N2,Ns) 



1 



— r- 

N 



2«iV n(a;) + l 



Figure 3: (Color online) This figure shows the amplitudes for the state, I^E'^), after three 
iterations of the algorithm. See the text for details. 



Substituting (j^ . this gives us a quadratic equation for B'': 

k k 

2-N{B^f + 2(^iV,/i,)s'= + Y,{N,- N,^^)Cl, -1 = 0, 
where we define A^^o = and 

k 

Cs,k = ^ ^ hj . 
j=s 

Solving this equation, and using the fact that > 0, we obtain 



(52) 



(53) 



2'^A^ 



2 

s,k 



2"A^ 



(54) 



This formula together with Eq. ()50p provide an explicit expression fl49|) for |\E''^) in terms 
of the numbers {A^^jJ^^ and {hs}s^i- 

To build the fcth feature, we apply the Grover operator G{Ok,tk) to the state We 
will now derive an expression for the integer "time" in terms of the features hi, . . . ,hk 
and the widths A^i, . . . , iV^. Given 1^^*^"^) as an initial state, let us define ^™ and 6™ to be 
the initial average amplitudes of the good and the bad states according to the oracle Ok- 



k-l 



(55) 



-^(aI-'Ni + At\N2 - iVi) + ■ ■ ■ + Atl{Nk-i - Nk-2) + B^-\Nk 



Nk-i) 



(56) 
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The initial separation, A™', between the good and the bad averages is therefore 

Af = g'r - W = ^^i^ . (57) 

Observe that developing a new feature of height hk is equivalent to increasing the initial 
separation A™ by hk. The final separation (after tk steps) between the good and bad 
averages is therefore 

Af = Af + h = ^^^^ • (58) 
Using (|37j) and (|3Hj) together with (jlSj) . we therefore have 



^ ( ■ ( 

tk = — arcsm 



A^ / Nk \ /Af /iVfc , , 

" ^ arcsin -^./t^ > (59) 



where 



and 



Uk = arccos ( 1 - ) , (60) 



2"A^ - A/'fc 

2 



iVfc(2»iV-iVfc) 



(61) 



To achieve a good fidelity between the state |\E'p) that we actually prepare and our "target" 
state |\l/p), we want the features hk to be as close as possible to the target values 6k 
defined in Eq. fl2Hj) . This motivates the formulas ()27flHn|l for tk which are obtained from 
the formulas (jSZHHH) by (i) replacing the features hk by the targets Sk, (ii) replacing the 
widths Nk by the measured values Nk, and (iii) by rounding to the nearest integer. 



4 Fidelity analysis 

In the above description of the algorithm, we have not specified how to choose the pa- 
rameters ?7c, rjg and a as a function of e (or, alternatively, of the initially given parameters 
A and t]). The optimal choice for these parameters depends on the estimates ni, . . . , ni/^ 
obtained in the quantum counting step. In this section we provide a rather generous worst 
case analysis which shows that the choice 

T], = eV54 , % = 0.99e2 , a = \log^ ^ - 3] (62) 

Vc 

guarantees, with probability greater than 1 — z/, the fidelity bound 

|(vI/,|vI/,)|>l-A. (63) 
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This bound is valid for arbitrary values of the n^. In most actual applications, much 
larger values of the accuracy parameters e, rig and r]c will be sufficient to guarantee this 
fidelity bound. 

We now show that Eq. (|3(jp. for the times which we motivated in the previous section, 
implies the fidelity bound (jUHj) under the assumption that the parameters rjc, rjg and a are 
chosen as in Eq. (jHH). 

Our starting point will be two sets of expressions for the tk, namely the definition of 
the tk, Eqs. (jTfHSUj) . in terms of the measured values Nk and the target values 6k, and 
Eqs. ()57flM|) above in terms of the actual values and h^. In subsection 14.11 we will 
derive an upper bound on the error \6k — hk\. This bound shows how accurate our 
algorithm is in achieving the target height, 6k, for the features hk- The overall accuracy 
of our algorithm, however, also depends on how accurate it is in achieving the correct 
width of the features. This accuracy is determined by the fraction of x values for which 
p'{x) 7^ p"{x) (see Figure Q). In subsection 14 . 21 we obtain an upper bound on this fraction. 
In subsection l4.31 we derive the fidelity bound and an upper bound on the probability 
that the algorithm fails due to a nonzero outcome of the measurement of the auxiliary 
qubits. And finally, in subsection 14.41 we show that the bound (fTT| on the phases (p{x) 
imphes the overall fidelity bound (fT^. 



4.1 Bound on \6k — hk\ 

It is convenient to perform the proof of the bound on \6k — hk\ in three steps. For this we 
note that tk depend on the values of 7™, and Uk- In subsection 14.1.11 we determine 
the range of possible values for 7™' and 'yf'^ that corresponds to the uncertainty in the 
measured values of Nk- Similarly in subsection 14 . 1 . 21 we determine the error range for ujk, 
and finally, in subsection 14. 1.31 we complete the proof of the bound. 



4.1.1 Error range for 7™ and 7^° 

Equation provides an explicit expression for tk in terms of {hk} and {Nk}: 

tk = F{{K},{Ns}). (64) 

In what follows it will be convenient to use auxiliary quantities {r^} defined as 

Tk = F{{6s},{Ns}). (65) 

The meaning of Tk becomes clear in comparison with (jMj) : Tk are the time intervals that 
correspond to the target features {6k}- Unlike {tk}, {Tk} are not necessarily integers. 

It will be convenient to rewrite the definition of tk given in Eqs. ()27riSm) in a slightly 
modified form: 

11 



tk 

where we use the following definitions: 



- + — (arcsm7^ — arcsm7^ 
2 iOk 



Uk = arccos ( 1 - ^ ) , (67) 
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Oik V S-^A^ afc V 2" AT 



and 



where 



k 

Bs,k = ^^j- 



It is also convenient to define 

fin Ei=l ^sSs , ini J2'lZlNJ, 

7fc = /^„.r.r and - 



In this notation Eq. ()65p can be rewritten as 



tUfcTfc = arcsin7^'^ — arcsin7™ 



Directly from the definitions we have 



it 



(Et; NsSsf + Nk{l- ElZliNs - Ns.,)Bl,_, 



By direct calculation we get 

fc-l fc-l k 



s=l s=l j=s-\-l 

k—1 k 



s=l " j=s+l 



The ± notation is an abbreviation for a double inequality (see the appendix) 

Ns > {rig — rjc)N we obtain 

fc-l k-l 

^ ilq ijc-' ^ 

s=l " s=l 

and similarly 
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Since t]c/rjg < 1/10, we therefore have 



(l±^Z^)3/2 



7r = 7r(l±10-). (79) 



and similarly 

7r = 7rfi±ic 

In the above formulas we have used the mean value theorem that states that for any 
function / that is continuous on the interval |x| < c, where c is some constant, and 
differentiable on |x| < c we can write 

/(x) = /(0)+x/'(xo), (80) 

where /'{xq) denotes the derivative of / at some point |a;o| < \x\. This gives, for example, 
that for |a;| < 1/10 

1 

1 ±2x, 



l±x 

(l±x)^/^ = l±1.6x. (81) 

The error bounds (fTSjl and (fTHj) were obtained by simplifying the somewhat tighter but 
unwieldy bounds using these methods. 



4.1.2 Error range for cuk 

The aim of this subsection is to determine the ratio between ujk and the true value Uk 
given by Equations (jHTjl and (|Hn|l respectively. Using the mean value theorem we have, 
by definition. 



uJk = arccos 1 — ( 1 ± 



%-Vc 



a.ccos|l^||U ' (±'^^U (82) 



1 - |1 - + 1,)!' ^ 2»iV„,-,,, 



where u is a real number such that \u\ < i]c/i]g. This implies 



- i|(i + «) - 



Since rjc/rjg < 1/10 we have 



^^<^, ^/r^>^^, ./2-^(l + u) > 1, (84) 
V9-Vc 0.%' ^ 0.9^2 V 2"iV^ ' ' ^ ' 



and therefore 



cDfc = cjfe±2^W^. (85) 
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It now remains to find a bound on J that is linear in ujk- We have, by definition, 



— — - = 1 — COS uJk = 2, sm — 
2^N 2 



Since > sin^ a; we obtain 



and therefore 



2''N - 2 



Vg 



4.1.3 Proof of the bound 

During the kth stage our algorithm creates a feature of height 

■ sm{uJktk - 6) - sin(-^fc 



hk = Oik 



where ^k is some initial phase. On the other hand, the target height 6k is 



Sk = Okd -j^ ysm{uJkTk - ^k) - sin(-^fc 



Hence 



\hi. - Sk 



Oik 



< 2au 



2«A^ 



2"A^ 
"iVT 



sin(t^fctfc - ik) - si\i{uJkTk - ik) 
^k{tk - Tk) 



sm 



/2«A^ I I 
< afeA/ t^fc^fc - Tfcl • 



Directly from the definition we obtain 

. f ■ -fin • ~im\ i 

ooktk = — arcsm — arcsm 7;. 1 ± cijfe 

UJk\ ' 



(86) 



^7) 



^9) 



(90) 



(91) 



(92) 



Using Eqs. (fTHj). (f?^ . and (jHHj) . using the fact that lOrjc/rig < 1/4 and using the inequal- 
ity ()151|) we have 

uJktk = . , \ fa±4,/l0^) ±u;fc. (93) 



(94) 



Now using 1)8111 and the fact that |Cfc| < tt we derive from the above equation 

UJkitk - Tk) = ±8— ± Uk . 

Vg 

Since a > 3 [see Eqs.flB^ and (0)], we have that 2Nk/ {2°'N) < 1/4, and we can therefore 
use the inequality p52|) to show that 



Uk = ±2 



2Nk 
2^N 



(95) 
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Since a < logg (r/^/r/c) — 3, we obtain the bound 

Mtk - Tk) = ±IQ- . (96) 

Vg 

The maximum possible value of the average amplitude of "bad" states is 

max 6 = , , (97) 

and the maximum possible value of the average amplitude of "good" states is 

max^ = (98) 
y/riN 



One can therefore write 



2 / / T^2 , (maX^)2A^; 

ai < (maxo) + 



2"iV - Nk 
1 1 

- 2^ ^ r/(2" - 1)A^ 

1 2 

< + 

< (w) 

Using this bound together with (j96|) we obtain for the error (|91|): 

32 n 

\hk - 6k\ < . ^ = ■ . (100) 

Substituting the parameters (jH^ into the right-hand side of the above inequality one can 
show that ^ 

\hk-6k\<^=. (101) 



4.2 Number of exceptional values 

When describing our algorithm in section |21 we have introduced functions p' and p" to 
distinguish two different approximations to the target function p. Namely, if p' is an 
approximation of p which is defined by the oracles Ok, then p" also takes into account the 
fact that we may not know the exact values of rifc. In our algorithm we therefore use p" as 
our target function which coincides with p' everywhere apart for a small fraction of values 
of X for which p'{x) ^ p"{x). In this section we obtain an upper bound on this fraction 
which will then be used in the next section where we derive bounds on the fidelity and 
the failure probability. 

For all j, we have \nj — nj\ < rjcN. For j < fi, we have fij < rjgN, and hence 

% < iVg + Vc)N . (102) 
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We now consider, for each k G {2, . . . , T}, all values of j such that < j < fk- For 
these values of j, we have 

hj < hf^_^ . (103) 

Since Uj > nf^,_^, we have 

1% - J < 2r/,iV , (104) 

hence 

\n,-Nk^i\<2r],N , (105) 

and finally 

\n, - Nk-i\ < Ar],N . (106) 
Since T < 1/e, we find that p"{x) ^ p'{x) for at most fiN values, where 



(107) 



4.3 Fidelity bound and the failure probability 

For k = T, Eq. ()49|) can be rewritten in the form 



x=0 j=l 



where 



Let us define 



Cj{x) 



1 if X < Nj 
otherwise 



Using this definition we have 



> 



J2 Vp(^[B^ + J2c,ix)h, 

x=Q j=l 

N-1 T 

5^V^(5^ + 5^c,(x)5^ 



N-l 



x=0 



(108) 



(109) 



(110) 



:iii^ 



where we have used the fact that p{x) = for x > N. Using pOlll . and since T < 1/e we 
obtain 



\d{x)\ < 



(112) 



Because ^p{x) < 1/ \/vN, this implies 



N-l T 
x=0 j=l 



€ 
V 



(113) 
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Rewriting Eq. (pUj) in terms of the coefficients Cj{x), 

T 



we can write 



N-l 



x=0 



In e 



(114) 



(115) 



For all X with a possible exception of at most nN values p"{x) = p'{x) (see Sec. I4.2|) . Let 
Se be the set of exceptional values of x for which p"{x) 7^ p'{x). We have 



N-l 



|(*p|*^)|>$^v^v^-Ev^v9R-v/pn^ -IS^IW . (116) 



x=0 



rj rj 



Since \/p{x), \/p'{x) and \/jf{x) are all bounded from above by 1/ ^JrjN we obtain 



> E V^X^- - ^ ■ (117) 



x=0 



rj rj 



By definition of p' 



(118) 



and since p is normalized we get 



y ?7 ?7 



(119) 



In order to continue we need to calculate \B'^\. This can be done by examining the 
normalization condition (^'^l^'-^) = 1. This condition reads 



2"Ar-i T 2 

[B^ + Y',{x)h) =1 



(120) 



x=0 j=l 

Using (jllUj) and (jll4|) this can be rewritten as 



or 



x=0 



Y + B"" + d{x)'^ -A = l 



(121) 



(122) 



x=0 



where 



A = J2 ([V¥{x)+B^ + d{x)y- (^^p"{x) + B^ + d{x)yj 

X&e ^ ^ 
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+2^ (123) 



Let us define 

e(x) = V^-V^- (124) 

Since p is normalized, equation p22|l gives a quadratic equation for B^: 



B'^f + 2t/5^ + 1/ = , (125) 



wliere 



and 



1 /^""^ \ 



1 / 

^ =^ E ( 2VpM(^(^) + e(x)) + (t/(x) + e(x))" 

+ ^E(P"(^)-P'(^) + 2^(^)(V^-V?M)) • (12T) 



Since \/p'{x) and ^/jfjx) are bounded from above by 1/ v^r/iV and since Se contains at 
most fiN elements (see Sec. I4.2p . we obtain with the help of Eq. ()112|) 

Similarly, since |e(x)| < e/y/rjN we have 

< ^i±i!!±if . (129) 
Since ^ < and 2" > see Eqs. ()107|) and (jHH) respectively, we obtain 

l^'^l < \u\ + ^/W+\v\ < 4i • (130) 

Together with Eq. ()119p this gives the lower bound on the fidelity, 

l(^pl^^)l>l-^, (131) 

where we have observed that e < 1/3 and used the bounds 2^* < /i < e^, which follow 
from our settings given in Eq. (jHH). The failure probability is 

Pfan = {2^N - 1)N\B^\' < ^ < lOA . (132) 
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4.4 Introduction of phases and the fidehty bound 

We now show that the choice A' = together with the inequahty (fTT|) . i.e., |0(x) 
< e'/^) imphes the overall fidelity bound (fT^ . The proof is straightforward. 

|(^r|^)| = [ y^p{x)p{x) exp[27ri(0(x) - 0(a:))] 

> A/p(a;)p(x) cos[0(a;) — (j){x)] 

X 

a; 

> 5^ (1 - e'VS) 
= lf*pl^p)l(l-AO 

> 1 - A - y . (133) 

5 Resources 

In this section we provide worst case upper bounds on the resources required by the 
algorithm. We distinguish between the resources that are needed for the state preparation 
part of the algorithm (subsection I5.1|l and the resources that are needed by the quantum 
counting that precedes the actual state preparation (subsection I5.2|l . 

5.1 Resources needed for state preparation 

5.1.1 Auxiliary qubits 

From our settings ()62j) we obtain 

^ < 54/6^ . (134) 

Vc 

We thus obtain for the number of auxiliary qubits 

a < \og^i%/Vc) - 3 < 3 + log2 e-^ . (135) 

5.1.2 Oracle calls 

Here we give an upper bound on the time resources needed by the algorithm. The con- 
struction of one feature requires at most 

max(tfc) < — (136) 
oracle calls. Using inequality (jHTj) we can therefore write 



maxtfc<vrY^^. (137) 
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From ()62|) we have 

V,-Vc> ■ (138) 

Since there are at most 1/e features and because A^^^ < {rjg — rjc)N and 2"- < we 
therefore have that the total number of oracle calls, noracie, satisfies the bound 

^Oracle < < • (139) 



5.2 Resources needed for counting 
5.2.1 Counting accuracy 

Consider an oracle O on the set of possible values of x. Using standard techniques 
we can count the number M of solutions of O within the absolute error AM 

AM < (V¥nM + ^ ,, ')2""^ , (140) 

where m is the number of auxiliary qubits needed by the standard quantum counting 
routine [TH] . 

We want AM < rj^-N , where rjc is the counting accuracy introduced earlier. This connects 
the desired counting accuracy rj^ with the number of auxiliary qubits m, 

where A = 2"™. Solving this equation for A in the case A > 0, we have 

A = 21-^^/2 (v/^qrr^-^), (142) 

where y = M/N. We see that, the bigger the value of a, the bigger m has to be in order 
to give the required counting accuracy r/c. We therefore set a to the minimum, i.e., a = 1 
(this doubles the range of x values to ensure reliable counting, see e.g. |T^). 

It is easy to check that the dependence of A on ?/ is monotonic. As we vary M in 
the range to — 1, the corresponding values of A vary between the limits y/2r]^ and 
\/2(\/l + Vc — 1) > Vcf^- It follows that the required number of auxiliary working qubits 
needed for counting with accuracy t]c is m < log?7~^. Thus we choose m = logr^"^. This 
choice guarantees the required accuracy of counting irrespective of the true value of M. 



(141) 



5.2.2 Counting probability 

The above counting procedure does not output the correct result with probability 1. For 
the procedure to work correctly with probability 1 — u we have to increase the number of 
auxiliary qubits from m to Oc which is given by 

a, = m + log2(2 + ^) = log2 • (143) 
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The number, Ai'count, of oracle calls that is required by the counting procedure is 



iVcount = 2"^ - 1 . (144) 

Substituting m = log2?7~^ and using Eq. (jU^ we obtain 

27(1 + 4z/) , , 

a, < log2 ^ , ^ . 145 

Since there are at most 1/e features the total number of oracle calls needed by the counting 
stage of our algorithm is bounded as 

iV'-l<^<?^. (146) 



6 Summary and conclusions 

In conclusion, we have described a quantum algorithm to prepare an arbitrary state of 
a quantum register of log2 N qubits, provided the state is initially given in the form of 
a classical algorithm to compute the complex amplitudes defining the state. For an 
important class of states, the algorithm is efficient in the sense of requiring numbers of 
oracle calls and additional gate operations that are polynomial in the number of qubits. 
The following table lists, for each stage of the algorithm, upper bounds on the number of 
oracle calls and the number of auxiliary qubits needed. 

oracle calls auxiliary qubits 

i.- 27{l+4i/) , 27(l+4i/) 

countmg loga Ve^ 

preparing ^ 3 + 3 logg ^ 

preparing |\E^) p- 

The bounds are not tight and can be improved by a more detailed error analysis. The 
total number of quantum gate operations depends on the implementation of the oracles. 
It is proportional to the number of oracle calls times a factor polynomial in log2 N if the 
functions p{x) and 4>{x) can be efficiently computed classically. 

Depending on the nature of the function p{x) and the prior information about p{x), the 
algorithm we have described in this paper can be optimized in a number of ways. For 
instance, the counting stage is the most expensive in terms of both oracle calls and ad- 
ditional qubits. If for some reason the numbers characterizing the oracles are known 
in advance, the counting stage can be omitted, leading to considerable savings. Further- 
more, in this case the fidelity bound can be guaranteed with probability 1, i.e., we can 
set 1/ = 0. 

In some cases the algorithm can be simplified if, instead of using the oracles defined in 
Eq. (fT^ . one uses oracles that return the fc-th bit of the expression p[x) / rjN . The 
general conclusions of the paper continue to hold for this variant of the algorithm, which 
we analyze in detail in Ref. pUj . 

Finally, by using generalizations of Grover's algorithm in which the oracles and the in- 
version about the mean introduce complex phase factors f7| ^] it is possible to reduce 
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the number of auxiliary qubits needed in the preparation stage of the algorithm. This 
leads to a reduction in the number of required oracle calls, and could also be important 
in implementations where the number of qubits is the main limiting factor. 
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A Appendix 

A.l Notation for double inequalities 

In this paper we have made a frequent use of the following convention. Let a, b and c be 
three numbers. The notation 

a = 6 ± c (147) 
is then understood to be equivalent to the double inequality 

b-c<a<b + c . (148) 

Furthermore, let g, e and F be functions. The notation 

hix) = J2F{9{x)±e{x)) (149) 

is then equivalent to the statement that h{x) can be written in the form 

Mx) = $^F(/(x)) , (150) 

where f{x) = g{x) ± e{x) for all x El. 

A. 2 Trigonometric inequalities 

Here we prove the following inequalities 

I arcsin(a; + t/) — arcsin(a;)| < 2Ay]t^, |i/|<l/4, (151) 

and 

I arccos(a: + z/) — arccos(a;)| < 2a/z/ , |t^|<l/4. (152) 
Consider the case > 0, which implies that for any x 

arcsin(x + u) — arcsin(x) > . (153) 

By inspection of arcsin- function we have that for < i/ < l/A the maximum value of the 
difference (arcsin(a; + i/) — arcsin(x)) is achieved for x — —1: 

max ^arcsin(a; + u) — arcsin(a;)^ = arcsin(i/ — 1) — arcsin(— 1) . (154) 
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In the case oi y < z the following equality holds 

arcsin z — arcsin y = arccos ^a/1 — y'^Vl — z"^ + yz^ . (155) 

Applying this equality to the right hand side of p54j) we obtain 

arcsin(a; + z/) — arcsin(x) < arccos(l — u) . (156) 
Let us now look for a constant c such that 

arccos(l — z^) < c^/h' . (157) 
Since arccos is a decreasing function the above requirement is equivalent to 

1 - z/ > cos(cv^) . (158) 
According to the mean value theorem, there exists u < Ca/z/ such that 

cos(cv^) = 1 ^cos(-u) , (159) 

and therefore the requirement ()158|) can be rewritten as 

1 < cos(m) . (160) 

It is clear that this requirement is guaranteed to be satisfied if we set c = 2. Indeed, 
2cosu > 1 for any nonnegative u < 7r/3 which includes all possible values of u that can 
correspond to c = 2 and v < 1/4. Since c = 2 guarantees that (jl57p is satisfied, we obtain 
from (|156p 

arcsin(x + z/) — arcsin(x) < 1\fv ■, z/ > 0. (161) 
The case of negative v can be treated in an analogous fashion leading to the inequality 

arcsin(x) — arcsin(x + z/) < 'l\/\u\^ z/ < 0. (162) 

The required inequality ()151|) follows trivially. Moreover, since 

arcsin X + arccos X = 7r/2, (163) 

we also obtain p52|) as required. 
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